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DECAY OF SOLUTIONS TO THE MAXWELL EQUATION ON THE 
SCHWARZSCHILD BACKGROUND 


LARS ANDERSSON, THOMAS BACKDAHL, AND PIETER BLUE 

Abstract. A new Morawetz or integrated local energy decay estimate for Maxwell test fields 
on the exterior of a Schwarzschild black hole spacetime is proved. The proof makes use of a 
new superenergy tensor defined in terms of the Maxwell field and its first derivatives. The 
superenergy tensor, although not conserved, yields a conserved higher order energy current 
The tensor vanishes tor the static Coulomb field, and the Morawetz estimate 
proved here therefore yields integrated decay for the Maxwell field to the Coulomb solution on 
the Schwarzschild exterior. 


1. Introduction 

The exterior region of the Schwarzschild spacetime is given in Schwarzschild coordinates by 
{x°‘) = {t, r,0,(j)) G R X (2M, oo) x [0, tt] x [0, 27r) with metric 

gabdx°'dx^ = (1 ~ 2M/r)dt^ ~ (1 ~ 2M/r)~^dr^ — r^{d9^ + sin^ Odtj)^). 

The exterior region can be extended, but for simplicity we will not treat the extension here. 
The Schwarzschild metric gab is static and spherically symmetric. In particular, the vector field 
= (5t)“ is Killing, timelike in the exterior, and orthogonal to the level sets of t. With our 
choice of signatureo the timelike condition is that > 0. In fact, as r —>• oo, we find t 1- 

Recall the condition for to be a Killing vector is ^(a^b) = 0. A related condition for a 2-form 
Yab to be a Killing-Yano tensor is ^(aYb)c = 0. The Schwarzschild exterior admits a Killing-Yano 
tensor given by 

Yab = -2r^smedeiad(j)b]. ( 1 - 1 ) 

A Maxwell test field is a real 2-form Fab satisfying the equations V^Fab = 0,\7[aFbc] = 0. From 
now on, unless otherwise stated we shall assume that Fab is a Maxwell field. 

To state our main result, we introduce the real I-form 

Ua = -r“^Var, 

and define in terms of Ua, Fab and Yab, a real 2-form Zab, and a complex I-form Pa by 


Zab = - ( 1 . 2 ) 

Pa= - \U^Zab - \iU^*Zab " - \lVb*Za\ (1.3) 

The superenergy tensors (in the terminology of [53]) for Pa and Zab are 

Haf, = P{ah) - yabP^Pc, (1-4) 

W,f, = - \Za^Zbc + IgabZadZ^'^. (1.5) 


Remark 1. The tensors Hob, Wat are quadratic in Fab cmd its derivatives up to first order (zeroth 
order for 'Wab)- Recall that a tensor Sab satisfies the dominant energy condition if SabV°‘Cf’ > 0 
for future-directed timelike vectors (^“. It follows from their definition that both H^b and Wab 
satisfy the dominant energy condition, cf. jT] 124] . The notion of superenergy tensor does not 
necessarily imply any conservation property. 

Given a vector field the superenergy current defined by H^b with respect to is H^bU^. 
For the static Killing field fP we have that the current H^b^^ is conserved, 

V“ (H.b^") = 0, (1.6) 

^In this paper, we follow the conventions of m and assume all fields are smooth unless otherwise stated. 
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see Lemma [TTl In contrast to the standard symmetric energy-momentum tensor for the Maxwell 
field, the tensors Hab and Wab are not conserved in general, ^ 0, V“Waf, ^ 0. On 

spacetimes admitting a Killing spinor satisfying an aligned matter condition, we have introduced 
another conserved tensor [T]. 

For a vector field z/“, a spacelike hypersurface S with future-directed normal iV“, and induced 
volume element /te, we define 

F;,(S) = (1.7) 

It follows from the dominant energy condition that if is timelike and future-directed, then 
Ev{'£‘) > 0. By the conservation identity (II. 61) . if Si, S 2 are hypersurfaces bounding a spacetime 
region O, then the superenergy with respect to on Ei and E 2 are equal. 

Let r“, be the orthonormal frame 

f“ = (I - 2M/r)-^/^{dtT, = r~^{deT, 

= r-icsc6»(a0)“, = (I - 2M/ry/^{drT, 

adapted to the foliation by level sets of t. Then if S* is a level set of t, the superenergy with 
respect to on St is given by 

= (|/3rP +1/3x1'+ l/3pP + l%n^'sin0drd0d0, (1.8) 

J 

where for a vector v°‘, = l3aV°‘- Similarly, 

Wff = id^fxl' + + \Zyz\^)- 

Hence, it controls all components of Zab, because the structure of Yab gives Zf-^ = 0, Z^^ = 0, 
see discussion below. 

We are now ready to state our main result. 

Theorem 2 (Energy bound and Morawetz estimate). Let Si and S 2 be spherically symmetric 
spacelike hypersurfaces in the exterior region of the Schwarzschild spacetime such that S 2 lies in 
the future o/Si and S 2 U —Si is the oriented boundary of a spacetime region II. 

If the real 2-form Fab = E[ab] *5 a solution of the Maxwell equations 

V^Fab=0, V[aFab]=0, (1-9) 

on the Schwarzschild exterior, and Zab and Pa are defined by equations dUD-dESD, then 

£;^(S2) = Ej(Si), (1.10) 

C 9M 79 

\^<^\i-des+^\Zab\ldpn < yi3j(Si), (1.11) 

where F^(Yii) is the flux (11.71) associated with and \Pa\i,deg and \Zab \2 are, respectively, the 
degenerate 1-form norm of Pa and the 2-form norm of Zab defined by 




I l^deg 

\Zab\l = 




M{r-ZMf{r-2M) ^ 


■|/3^ 


Tl ’ 


2 _ (r - 2M) 


W 


TT- 


Remark 3. The approach presented here can be extended to the inhomogenous Maxwell equations. 


Estimates of the form (11.111) are called Morawetz or integrated local energy decay estimates, 
since they show that local norms in space are integrable in time. Since the right-hand side 
of inequality (ll.ll|) depends only on initial data, it follows from the Morawetz estimate that the 
energy in stationary regions must tend sequentially to zero for large times, which demonstrates 
that the held disperses. In the present discussion the energy refers to the superenergies dehned in 
terms of Pa and Zab^ and dispersion refers to the fact that these helds tend to zero asymptotically 
in any stationary region. Although the parts of the Maxwell held controlled by Zab disperse, the 
non-radiating Coulomb solution, which is a bound state, remains. Below, we discuss that the 
helds Pa and Zab vanish when evaluated on the Coulomb solution. 
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Recall that the standard Maxwell energy-momentum tensor, given up to a constant factor by 

Tab = -^Fa^Fbc + IgabFcdF^'^, 


is traceless and conserved if Fab is a Maxwell field, see [52J Chapter 5]. We remark that the 
analogue of Theorem [2] cannot hold if we replace Hab by Tab, due to the fact that the energy 
density defined in terms of Tab is non-zero when evaluated on the Coulomb solution. In fact, for 
a radial Morawetz vector held, it is not difficult to see that the bulk term V°‘{TabA^) = Tab'V°'A^ 
must change sign at r = 3M, see [31 section 7]. 

The main motivation for this work is to develop techniques that might be useful in approaching 
the Kerr stability conjecture. The Kerr family of spacetimes is a two-parameter family of solutions 
of the Einstein equation describing rotating black holes. The Schwarzschild spacetimes make up 
a one-parameter subfamily corresponding to zero rotation speed. The Kerr stability conjecture is 
that the Kerr family of solutions is asymptotically stable under the Einstein equation, although 
the Schwarzschild subfamily in isolation is not expected to be stable, except when restricting to 
the axially symmetric case with zero angular momentum. A natural strategy is to understand 
solutions of first the null geodesic equation, the wave equation. Maxwell equation, and then some 
reasonable linearisation of the Einstein equation, before approaching the full nonlinear Einstein 
equation. 

To compare Theorem|3]and its proof with earlier work in the literature, one must consider the 
null or spinor decomposition of a 2-form. To do so, one first introduces a complex null tetrad 
(Z“,n“,m“,m“) normalized so that gabFn}’ = —gab'm'^fh}’ = 1 and all other inner products are 
zero. For example, in the exterior of the Schwarzschild spacetime, one can use 


r = z“), 

V2 


n“ = ^(T“ - Z“), 
V2 


±(X^ + ^Y^). 


(1.13) 


This is a principal null tetrad. For a 2-form Fab, the Newman-Penrose scalars with respect to a 
null tetrad (not necessarily principal) are 


00 = FabFm^, 01 = ^{FabFn’’ + Fabm‘^m^), 02 = FabnFn^. 

The components 0o, 02 are called the extreme components and 0i is the middle component. The 
extreme components depend on the scaling of the tetrad, while the middle component is a true 
scalar. Since the vector fields dg and cscOd^ are not continuous at 0 € {0,7r}, the extreme 
components are not smooth as functions even if Fab is smooth, but are smooth when viewed as 
sections of an appropriate complex line bundle, cf. [22], see also [8]. 

A key tool in most previous work on this topic has been the existence of second-order, decoupled 
equations for individual null components. In the exterior of the Schwarzschild spacetime, the 
middle component 0i of the Maxwell field satisfies (after rescaling by a power of r) a wave 
equation with real potential, which is the Schwarzschild case of the Fackerell-Ipser equation 
[16] . In the Kerr spacetime, the potential in the Fackerell-Ipser equation is complex, and the 
extreme components satisfy a more complicated, complex second-order hyperbolic PDE known 
as the Teukolsky equation m- Both the Fackerell-Ipser and Teukolsky equations have analogues 
arising in certain linearisations of the Einstein equation, making them particularly interesting as 
part of the strategy for approaching the Kerr stability conjecture. In fact, the potential in the 
Fackerell-Ipser wave equation on the exterior of the Schwarzschild spacetime is a multiple of the 
potential in the Regge-Wheeler equation governing odd parity gravitational perturbations of the 
Schwarzschild metric, cf. ESjug. 

The Maxwell equation in the exterior of the Schwarzschild spacetime has stationary solutions 
given by 

Fab = ^{qET[aZb] + 


For these, the Maxwell scalars are respectively 0o = 0, 0i = {qe + iqB)l'r^, and 02 = 0. These 
are called the Coulomb solutions and are the only solutions that do not decay to zero. Thus, one 
commonly says the Coulomb solutions are supported in the middle component. Clearly, these 
must be excluded for a decay estimate to hold. 
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At each point, the middle component of the 2-form Zab is zero and its extreme components 
(modulo a sign) are r times the corresponding extreme component of the complex self-dual 2- 
form Fab + ii*F)ab- Thus, by working with Zab, we have introduced a variable that geometrically 
excludes the Coulomb solutions. Since at any point the components of Fab can be freely specified 
while constraining Fab to be smooth and compactly supported on any spacelike hypersurface Si, 
the estimate in Theorem [2] is non-vacuous. 

Decay for the Maxwell equation in the exterior region of a Schwarzschild spacetime or its 
generalisations has been studied in [HI 1 HU Ha 1 m]. It has long been clear that, since 
is a time-like Killing vector, there is a conserved positive energy for Fab, and that this can be 
used as a foundation for subsequent decay estimates. All works have used both a second-order 
PDE for certain components and the fact that such solutions arise from the original first-order 
Maxwell equation. The paper m is the only previous work we know of to prove decay using the 
Teukolsky equation; the conclusion of that work is that at fixed r the extreme components go to 
zero, and the method relies on an integral representatiorQ. In contrast, the remaining work has 
focused on the middle component, which satisfies the Fackerell-Ipser equation and which requires 
the contribution from the Coulomb solution to be explicitly subtracted off. 

Since the Fackerell-Ipser equation is a wave-like equation, we recall that the last 15 years have 
seen the proof of Morawetz and related decay estimates for the wave equation in the exterior of a 
Schwarzschild black hole [iniiniiiiiiii] essentially using vector-field methods and following ideas 
introduced for a model problem in m- In the Kerr case, because of the ergoregion generated by 
rotation, energy and Morawetz estimates had to be proved simultaneously instead of sequentially. 
This was first done in the very slowly rotating case [ia[26l[5] and more recently the full range 
of subextremal Kerr black holes M- Because of the complicated nature of the orbiting null 
geodesics around Kerr black holes, it was necessary to extend the vector-field method. In [HE] 
and |26j , this was done by blending it with separation-of-variable and operator-theory techniques 
respectively; whereas in [5] , the first and third authors were able to work with purely differential 
operators. This involved using hidden symmetries, second-order differential operators that take 
solutions to solutions but which are not decomposable into first-order symmetries. In [5], we have 
explored the analogue for the Maxwell equation. 

Returning to the Maxwell problem, |5] proves a Morawetz estimate for solutions to the Fackerell- 
Ipser equation on the exterior Schwarzschild spacetime and uses this to prove a Morawetz estimate 
and decay for the components of Fab at fixed r and a hierarchy of decay rates along null in¬ 
finity (where t —>■ oo but t — r remains bounded), and relies upon energies arising from the 
vector-field method. The paper m proves t ^ decay at fixed r for the middle component by 
summing Sigal-Soffer propagation estimates for separated modes, and |5S] treats inhomogeneous 
Maxwell equations on a very general class of stationary, spherically symmetric black-hole space- 
times and the results extend beyond the exterior region. The proof uses energies and vector-field 
methods first to prove a Morawetz estimate for the middle component and then to extend this 
to a Morawetz estimate for all components. The paper m proves a “black-box” result, in which 
a Morawetz estimate of the strong form proved in [25] is assumed and shown to be sufficient to 
imply a decay rate at fixed r without using the details of the proof of the Morawetz estimate. 
The paper m similarly proves a decay result based on the assumption of a Morawetz estimate. 

In [4], the first and third authors proved a Morawetz estimate outside a very slowly rotating 
(and hence non-spherically symmetric) Kerr black hole. So far, it has not been possible to prove 
decay estimates for the Maxwell equation without passing to a second-order equation. Energies 
for second-order equations are naturally at the F[^ level rather than the level for first-order 
equations. In the Schwarzschild case, the Fackerell-Ipser equation has a conserved, non-negative 
energy, so there was no obstacle. However, for the Fackerell-Ipser equation on the Kerr spacetime, 
the potential is complex, which destroys the conservation structure. This could be compensated 
for, but only by introducing fractional derivative operators (more precisely, fractional powers of 
the separation constants after separation of variables). It seems reasonable to imagine that the 


^We have been unable to verify the claims in dzi concerning the spin-2 case. In particular we have been unable 
to verify that the energy presented in |17l Eq. (1.4)], also discussed in appendix A of that paper, is conserved. For 
the Maxwell case, the analogous claim follows from the conserved positive energy for Fab- 
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original type energy could equally well be combined with a fractional derivative operator to 
prove a similar result. 

In this paper, we develop a new method for treating the Maxwell equation which does not 
require introducing a second-order equation. This opens the possibility of applying the differential 
hidden symmetry operators from [5] , which might allow for the proof of decay estimates in the Kerr 
case without using fractional derivative operators or fractional powers of separation constants. 

The principal variables in this paper are Zab and /3a, which are constructed from the extreme 
components, which satisfy the Teukolsky equations, rather than the middle component, which 
in the case we consider satisfies the Schwarzschild case of the Fackerell-Ipser wave equation. We 
believe this has two major advantages. First, it provides a local way to exclude the Coulomb 
solutions, without integrating over spheres. Second, we have recently found a new conservation 
law for the Maxwell equation [T] which generates level energies. Curiously, although Zab in 
this paper and in [T] are the same, the auxiliary 1-forms /3a and rja and the associated symmetric 
2-tensors Hat, and Vat, differ. Nonetheless, in the sense given in [T], the leading-order parts of 
Hat, and the conserved Vat, are the same. We are currently investigating how to replace Hat, by 
Vat, so that the argument in this paper can be extended to the rotating Kerr case. 

By applying the Maxwell equations, the variable /3a can be transformed from an expression 
in the derivatives of the extreme components to an expression purely in terms of the middle 
component. Recall that in CHP notation, the Maxwell field. Fab, is decomposed into (j)Q, (j)i, (/) 2 , 
and we refer to (j)i as the middle component. In particular, j3a = VaT — CaT, where T is r times 
the middle component of Fab- Similarly, one might view Hab as a type of energy-momentum tensor 
(albeit not conserved) for the Fackerell-Ipser equation. Following this, one might be tempted to 
conclude that Theorem [2] follows immediately from the Maxwell equation and the results about 
the middle component in [5] or [5S]. We emphasise that such a naive interpretation is false. First, 
Theorem [2] includes undifferentiated factors of 0o and 4)2, which arise from the components of Zab 
and which cannot be recovered from the middle component using local operators. Second, in [5], 
a Hardy estimate in the radial derivative was applied to the middle component; in Section [2.21 a 
Hardy estimate in angular derivatives is applied to the extreme components. Morawetz estimates 
typically require a Hardy estimate to obtain a globally positive coefficient for the undifferentiated 
terms. 

The argument in this paper very closely follows the argument used for the wave equation. 
From a vector field A°' and a scalar field q, we define the a current Pa in terms of Zab and Pa bjl 


Pa = HabH" + \qP^Zab + \qP^Zab - \iqP^*Zab + \iqP^*Zab 
+ ^ZbcZ^‘^Waq-kZa‘^ZbaV'’q. 


We now choose 


(r-3M)(r-2M) „ 

^ ^ WrJ , 

9M2(r-2M)(2r-3M) 

4j,5 


9 = 


(1.14) 


(1.15a) 

(1.15b) 


The motivation for this choice is explained in the proof of Theorem [THl The vector field is the 
same as occurs in the analysis of the wave equation. In particular, it has the geometric property 
that it points away from the orbiting null geodesics, which in the Schwarzschild spacetime are 
located at r = 3M. 

For any spacelike hypersurface E, we define 


[ HabC'-VM/rE, 
JY. 


(1.16) 


E^+A,qiZ.)= f (Hab^'+Pa) V“d/rE. 

Jy 


(1.17) 


^This can be compared with the corresponding expression for the wave equation, 


Pa = T[u]abA^ + quVaU - {Vaq)u‘^/2, 

where T[u]ab is the energy-momentum tensor for the wave equation. Thus, we use Hab place of the wave 
equation’s energy-momentum tensor, and we use the terms involving q in lITTill in place of qu\/aU — (Vaq)'i^^/2. 
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As discussed above, E^{T,) is nonnegative and conserved. In Section|4]we show that on spacelike, 
spherically symmetric hypersurfaces, the energies E^{'E) and £’j+a,(j(S) are uniformly equivalent. 
In Section [SJ we show that the integral over any spherically symmetric region of —VoP“ 
dominates the integral on the left hand side of (11.1111 . The proof of this fact relies on a spherical 
Hardy estimate, see Lemma [THl Finally these facts are combined in Section [5] to yield a proof of 
Theorem [2l 


2. Preliminaries 

For the remainder of this paper, we will make use of the 2-spinor formalism, following the 
conventions of |22) . The exterior Schwarzschild spacetime is oriented and globally hyperbolic and 
hence also spin. 

The spin group is SL(2, C) which has the inequivalent spinor representations and (C^. Un¬ 
primed upper case latin indices and their primed versions are used for sections of the corresponding 
spinor bundles, respectively. The correspondence between spinors and tensors makes it possible 
to translate all tensor expressions to spinor form. The action of SL(2, C) on leaves invariant 
the spin metric cab = £[ab] > which is used to raise and lower indices on tensors. The metric gab 
is related to cab by gab = ^ab^A'B', and the scalar and Weyl curvatures are represented by the 
spinor fields A, 'i’ abcd, respectively. We denote the space of symmetric spinors with k unprimed 
indices and I primed indices by 

The principal null tetrad E, n°^ ^ m°'^ fh°' corresponds to a principal spin dyad OA,bA via la = 
OAOA',na = LAlA,'>Tia = OA^A'jEia = laoa'■ Given a symmetric spinor field (pa---da'-- d'^ we de¬ 
note the dyad components by (pui where i,i' denote the number of contractions with , 

respectively. For example, a 1-form (fAA', has components ipoo' = '^AA’O'^d^ ,■ ■ ■ = 

A -A' 

ipAA'l ■ 

The 2-spinor formalism makes it possible to systematically exploit decompositions in terms 
of irreducible representations of the spin group SL(2,C). For the spinor calculations, we have 
made use of the SymManipulator package [5], developed by one of the authors (T.B.) for the 
Mathematica based symbolic differential geometry suite xAct |20) . 


2.1. Fundamental operators. 





^k,i 8>k-i,i+i o,nd SAk. 

.1 '■ Sk,l 

{^ky^)Ai.. 

A' 


_ \-7BB' A\ _ 


i'^k,l^)Ai. 

A' 


_ T-7 b'. A\ .. 

= V(^j ^PA2...Ak+i) 

B‘ 


A'. 




{%,1<P)Ai. 

A'. 

1 





The operator Sik,i only makes sense when k > 1 and / > 1. Likewise is defined only 
if 1 > 1 and only if fc > 1. From the definition it is clear that the complex conjugates of 

and are and {^i,k(p) respectively, with 

the appropriate indices. With these definitions, a Killing spinor of valence (/c, 1) is an element 
KA-F^'"’^' G her 

The fundamental operators appear naturally in the irreducible decomposition of the covariant 
derivative of a symmetric spinor field, see [U Lemma 15]. In most calculations we freely make 
use of such decompositions. We shall make use of the following commutator relations for the 
fundamental operators. 
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Lemma 5 ([2j Lemma 18]). Let ipAB S § 2 , 0 ; a scalar. The operators and 


satisfies the following commutator relations 

fi^iA%,o^)AB = 0, (2.1a) 

(f^l,%,o^)A'B'=0, (2.1b) 

('^3,i-^,o</’)abcd = ‘2''^{abc^Td)Fi (2-lc) 

i^3,1^2,0T)AB = — AB — SiAifAB + ABCDt'"^ ■ (2-ld) 


2.2. Hardy estimate. Recall the scalar components of spinors transform under tetrad rescalings 
as sections of certain complex line bundles. The GHP operators ]3,]3^,S, are defined on such 
scalars as the appropriate covariant derivative along the null tetrad legs Z“,n“, m“, respec¬ 
tively, see m for details. 


Lemma 6. On any sphere with constant r in the Schwarzschild spacetime, with tpo and ip 2 the 
extreme components of a smooth symmetric spinor field ipAB we have 



9 ' (fol'^dfiSr-, 


(2.2a) 





< 



9(/32pdMS;.- 


(2.2b) 


Proof. Expand ipo and ip 2 in terms of spin-weighted spherical harmonics (see |221 Section 4.15]) 

OO I 

E E iYi,m{0,4>), (2.3a) 

m— — l 

OO I 

T2{6,4>) = EE hl,m-lYl,m{0,4>). (2.3b) 

/ = 1 m——l 

From Eq. (4.15.106)] we have 

S'^o{S,4')= E E —oYi,m{d,(l)), (2.4a) 

S^2{e,fi)= Y. (2.4b) 

1—1 m— — l ^ ^ 

Through the orthogonality conditions [HI Eq. (4.15.99)] we get 


p OO 1 

/ l^opd/rs. = 47ry^ |az,mP, 

(2.5a) 

1—1 m——l 


p OO 1 

/ \(p2\^dpSr = E 

(2.5b) 

m=-l 


’St- ? 1 7 

^ 1—1 m——l 

(2.5c) 

1=1 m=-l 

(2.5d) 


This proves the desired inequalities. □ 


3. A FIRST ORDER SUPERENERGY FOR THE MAXWELL FIELD 

3.1. General spacetime. In this subsection we shall work in a general spacetime with a Killing 
spinor kab of valence (2,0) satisfying {^ 2 ,o>^)abca' = 0 and such that kab^^^ 0. Recall that 
a Maxwell field Fab corresponds to a symmetric spinor of valence (2,0) by 

Fab = eA'B'^AB + £Ab 4 >A'B'- 
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The source-free Maxwell equation takes the form {‘^2o4')aa' = 0. Instead of constructing energy- 
momentum tensors from the Maxwell field itself, we will use a first order expression. In the paper 
[T] we have introduced a conserved tensor for the Maxwell field on any spacetime admitting 
a valence (2,0) Killing spinor with aligned matter. Here we shall take a different approach 
and construct a superenergy tensor satisfying the dominant energy condition which is however 


not conserved. This tensor is constructed from a quantity Paac which can be defined for any 
spacetime with an algebraically general valence (2,0) Killing spinor. We remark that the leading 
order term in the tensor constructed here agrees with the tensor constructed in [T]. 

As we shall discuss later, the spinor fields introduced in the following definition correspond to 
the tensor fields given in the introduction. 

Definition 7. Assume that kab is a valence (2,0) Killing spinor, such that 
a solution to the source-free Maxwell equation ['^ 2 oP)aa' = 0, and define 

kcdfP^ 0, Pab 

O 

II 


(3.1a) 

UaA' = - ^^AAAog{-KCDK^^), 


(3.1b) 

T = k^^Pab, 


(3.1c) 

©AB = — 2K(a‘^</>B)C) 


(3.1d) 

PaA' = ('^ 2 ^, 0 ©)^^' 4 ” A'QaB- 


(3.1e) 

Lemma 8. Definition^ implies the relations 



Paa' = ~ Uaa'A^ -\- (■PopT)aa' 1 


(3.2a) 

(^ 1 , 1 / 3 ) = -U^^'Paac 


(3.2b) 

P^1.iP)aB = U(^A^ Pb)ac 


(3.2c) 

P^1aP)a'B' = U^(A'P\A\B')- 


(3.2d) 


Remark 9. Observe that we have not assumed anything about the matter or the behaviour of 
f,AA' ■ 


Proof. Expanding the definition of Uaa' yield 


Uaa' 


3{kcdk^^)' 


(3.3) 


We can write the Maxwell field in terms of Qab and T through the relation 

4>AB = — 'dA^KBC + Kab'^- (3.4) 

Differentiating the definition of T, using the Killing spinor equation and the Maxwell equation 
gives 


{^2fi(t>)ABCA' = A'<fAB + (=%,oT)aA'- (3.5) 

Expanding the definition of I3aa' , using (USD and (EH) we get 

PaA' = U^A'QaB — A'pAB + {^^2,Op) ABC A' 

= A'QaB — A'pAB + i%,o'^)AA' 

= -Uaa'T + {%^oVaa'- (3.6) 

Using the commutators (I2.1dll . (I2.1cll and the Killing spinor equation, we get 

(a^IsOaB = PiPi,i'^2.0^)aB = —QAkab + ABCDfi^^ , (3.7) 

0 = -^2,qk) ABCD = '^{ABC^^D)F- (3-8) 


Differentiating (1331) . using (13.4L the commutator (I2.1dll . the Maxwell equation, (13.81) and (13.7|) . 
we get a version of the Fackerell-Ipser equation 


(^i.i^o.oT) 


- 4AT + ^ABCD^i^^P^^ 


{’j^iApAB J. 

3{kcdk^°) 


(3.9) 
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Direct calculations using the Killing spinor equation yield 

(^i,i£/) = - Uaa'U^^' + ^ C^i.iC/)ab = 0, {^I^U)a'B'=Q- (3.10) 

The equation (I3.2bl) follows from (I3.2al) . (13.101) and (13.91) . The equations (I3.2cl) and (I3.2dl) follow 


from (I3.2al) . (13.101) and the commutators (I2.1al) and (I2.1bl) . □ 

Lemma 10. The components of (pAB, ^ab ond Paa' in a principal dyad are related by 

©0 = — 2ki^0; ©1 = 0; ©2 = 2Ki(f>2, (3.11) 

Poo' = S^©0) /3oi' = t>^©0) Pw' = ~t>©2j Piv = —9 ©2. (3.12) 

The superenergy tensors for Paa' and ©^is are given by 

^ABA'B' = \PaB>Pa'B T \Pba'PB'a^ (3.13a) 

^ ABA'B' = QabQa'B'- (3.13b) 

The notion of superenergy tensor extends to spinors of arbitrary valence, see |24) . 


We close this section by giving the correspondence between the spinor fields kab,&ab and 
the tensor helds Yab, Zab- We have 

Yab = ^ii^A'B'KAB — ^ABKA'B'), 

Zab = ^A'B'QaB + ^AbQa'B'- 

The normalization of Yab is chosen for convenience. For the Schwarzschild spacetime with a 
principal dyad (oa, ^a) we have 

kab = lro(^A^B)- (3.14) 

In the Schwarzschild case, the above definitions with kab given by (13.141) yield the fields Yab, 
Ua, Zab, Pa, Haf,, Waf, as in the introduction. Finally, the Morawetz current Paa' given in 
tensor form by (11.141) can be written in spinor form as 

PaA' = HabA'B'A^^' - IqpA'^^AB - ^Pa^'Qa'B' + A^QA'^\%m)BB'■ (3.15) 


3.2. Schwarzschild spacetime. 


Lemma 11. For the Schwarzschild spacetime we have 

W^^'HabA'B' = -Uaa'P^^'Pb'b, 
= 0 . 


In particular, TiABA'B' is a positive conserved current. 


(3.16a) 

(3.16b) 


Proof. For the Schwarzschild spacetime with a principal dyad {oa,i-a) the Killing spinor kab is 
given by (I3.14|) . This gives immediately 


f,AA' 


Uaa' 




(r 


-r Vaa'T 


- 2M)i/2f 
H72 ’ 

(r - 2MY/^Zaa' 

j,3/2 


(3.17a) 

(3.17b) 


Computing the divergence of (j3.13aL doing an irreducible decomposition of the derivatives, 
using Lemma [8] and using the reality of Uaa' gives (I3.16a|) . Equation (I3.16bl) , then follows from 
i^^'UA'A = 0 . □ 


Remark 12. For the Kerr spacetime with non-vanishing angular momentum, the 1-form Uaa' 
fails to be real and the current is not eonserved. 
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4. Positive energy 


Before we do any integrated decay estimates, we will verify that the energy (11.171) will be 
positive and uniformly equivalent to the energy (11.161) . 

Theorem 13. Let and q be given by (ll.lSal) and (jl.lSbl) . For any constant |ci| < 10/9 

and any spherically symmetric slice S with future pointing timelike normal such that 

N^^'Naa' = 1 we have a positive energy 


+Ci-PAA')dp^^ > 0. 


(4.1) 

(4.2) 


Proof. We denote the superenergy tensor of a vector field by 

^F‘l\aAA']ABA' B' = \aAB'aA'B + ^OiBA'^B' A- 

This tensor satishes the dominant energy condition. 

Now, as the slice is spherically symmetric, the future pointing timelike vector field is 

spanned by and with coefficients depending only on r. We can in general write it as 

= {w{r) + \w{r)-^)f^^' + {w{r) - \w{r)-^)Z^^\ (4.3) 

for some radial function w{r) > 0. We can then verify the identity 

FN^^'CPa'^&ab + Pa^'Qa'b') = 2e-^llABA'B'N^^'f^^' + 2eN^^'f^^'Q abQa'B' 

- 2e-iiV^^'f^®'SEi[/3AA' T bf^eN^ a'QacUba’B' 

-bf^ef^^'f^^'QABeA'B', (4.4) 

where bi = Taa' = jwfr) ^ + w{r) > 0. The dominant energy condition then gives the 
Cauchy-Schwarz inequalities 

±N^^\Pa'^Qab + Pa^'Qa'B') < 2e-^llABA'B>N^^'f^^' + 2eN^^'f^^'QABQA'B'. (4.5) 
Using this in the definition (13.151) and the positivity of g, we get 
iiV^^'PAA' > liABA'B'N^'^’iFA^^' - e-^qf^^’) 

+ N^^'QABQA'B'{-eqf^^' ± \{%pq)^^'). (4.6) 

Expanding in the dyad gives 


N'^^'fBB'e^^e^,^, = u;(r)|0op + 

N^a'^bb'qabQa'B' = iT(r)|0oP - 


1021^ 

Fwir) 

102^ ^ 

Aw{r)' 


N^^'fBB'nABA'B' = iT(r)| g' 0o|" + u;(r)| ^ 0o|" + 


i>02p , |g02p 


Fwir) Fw(r) 


' |2 


> u;(r)| g 02 


|g02p 


Fwir) 

Hence, from (I2.2al) and (I2.2bl) we get the Hardy estimates 


N^^'fBB'Q^^QA’B'dps. < ^ 


|iV^^'Z^^'0AB0A'B'MMS. < y / N^^'T^^^'Haba'bApls^ 


N^^'fBB'UABA'BApiSr 




Hence, we can make the estimates 


where 


[ ±N^^'PAA'dpSr> [ liABA'B'N^^'{±ABB' 

JSr- JSr- 

^AA' ^ ^-iqfAA' ^ d\q'{r)\r^/^{r - 2Mf/‘^f^^'. 


(4.7a) 

(4.7b) 

(4.7c) 

(4.8a) 

(4.8b) 

(4.9) 

(4.10) 
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Here we have used the ansatz q = q{r), which yield (^,o'?)aa' = — (t -2MY^\-^^Wr)ZAA'- 
With our choice of q{r) and 

we get 

bAA' ^ |M2r-4(6M2-13Mr + 6 r 2 + |r-3M|(3r-5M))C^'^'. 

In the region r > 2M, we have 

|M2r-^(6M2-13Mr + 6r2 + |r-3M|(3r-5M)) < 

Hence 

/ VAA'^^J-Sr ^ [ '^aba'B'N^'^ {±A^^ — )dfiSr- 

JSr JSr 

This gives 

[ N^^'{UabA>B'^^^' +CiPAA')dflSr> [ H^BA'BW^^'((l-||ci|)e^^'+CiH^^')d/is.. 
JSr JSr 

(4.15) 

With our value of /(r), the vector field + C 2 A^^ is future pointing and timelike if |c 2 | < 2. 
Hence, the right hand side of (I4.15P is non-negative if |ci| < 10/9. □ 

Corollary 14. For any spherically symmetric slice E with future pointing timelike normal 
such that Naa' = 1 the energies E^(T,) and E^+A,q{E) are uniformly equivalent, 

^e^{j:) < E^+aA^) < (4-16) 


(4.11) 

(4.12) 

(4.13) 

(4.14) 


5. Integrated decay estimate 
In this section, we will prove an integrated decay estimate. 

Lemma 15. On any sphere with constant r, we have 

[ <^/ \/3f\^ + \^3d"d^iSr■ ( 5 . 1 ) 

JSr ^ JSr 

Proof. We have 

l/3zP + l/9zP = |302p + |S'0oP, (5.2a) 

Wff = T^^'t^^'QabQa'B' = ||0oP + 5|02p. (5.2b) 

The equations (I5.2ap , (I5.2bl) and Lemma IH gives the desired inequality. □ 

To analyse the positivity of the bulk term, we will write it in terms of quadratic forms. We 
shall use the following notation. 

Definition 16. For any vector field we define the quadratic form 

Ei[bo,bi,b2,BAA'] = {2b2f^A'f^B'F{bl-2b2)Z^A'Z^B')uU^'i'^'A) 


- {bo + jbi)L'AA'A''^, 

(5.3) 

where bo, bi, 62 are scalar fields. 


In particular. 


£’i[— 5 , 2, i, (3aa'] = /3f!P + /3f P- 

(5.4) 

Lemma 17. Ei[bo, bi, b 2 , baa'] > 0 for all if and only if 


0 < &i. 

(5.5a) 

0 < 62 , 

(5.5b) 

max(— 62 , —5i + ^ 2 ) <bo< 62 - 

(5.5c) 
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Proof. The eigenvalues of the quadratic form Ei[bo,bi,b 2 ,VAA'] at a point can be calculated by 
expanding in i^aa' in dyad components. One finds the eigenvalues 


^ 0+^1 — ^ 2 , —^ 0 +^ 2 , bo + b2, &0 + & 2 - 

Requiring each eigenvalue in (I5.6|) to be non-negative gives the result. 


(5.6) 

□ 


Theorem 18. Let Paa' be given by (13.1511 with the choices (I1.15a|) and (Il.l5bl) . Then we have 
the estimate 


r r \ ]\{ 

J^-{^lsP)dHQ > J^-\PAA'\l,deg+ 

for any spherically symmetric spacetime region O of the Schwarzschild spacetime. 


(5.7) 


Proof. In order to motivate the choices (I1.15all and (Il.l5b|) , we shall start by considering general 
radial and q. From the form (13.151) . the definition of /3aa' and the properties from Lemma [ 5 ] 
we get 

-(^l,lP) = - P^'^{^,aA)aba'B' + Pa'a{\{^iaA) + A^^'Ubb' - q) 


+ eABQA'BiU^^'{%.gq)^^' - ^') ■ 


BB' 


\ABA'B’\ 


With the ansatz 


^AA' _ \AA' _ firy^'^Z^^' 


A^^ =f{r)idrr^ = 


{r - 2 M)i /2 ’ 


q = q{r), 


we get 


^ + /'(r), 


(■^,i^)ab 


UB' ^ /(r)(-(r - + (r - Myy^'Zsf)) 

{r — 2M)r 

-Z^A^^'z^f^fir), 

(r - 2My^ZAA>q'ir) 


iyi,0<l)AA' — — 


r.1/2 


{^i,i-%,oq)ab 


UB' _ (r - 3M)(f(^(^ Tb)^'^ - Zy^'ZB)^'y'{r) 


+ 


(r - 2M)Zy^'ZBf'h"{r) 


(5.8) 

(5.9) 

(5.10a) 

(5.10b) 

(5.10c) 

(5.10d) 


We can now write the divergence of Jaa> in terms of the quadratic forms Ei and Wff, making 
use of the fact that since the middle component of Qab is zero, 'Wff = 'W gg. This gives 

‘^Mf{r) f{r)ir - 3M) 


/(r)(r - M) _ 
2r{r-2M) 


-(^i,iP) = ii;i[g(r) + - i/'(r ),+ f^r), ^ ^yy ^A' 


(r — 2M)r 

{r-2M)q'lr) [r-2M)q"{r). 

~) • 


2{r-2M)r 


(5.11) 


Later we will need the Hardy estimate in the form of Lemma[T2]to handle a negative contribution 
from the term. Therefore, we make use of (1^ and extract a term g{r) (|/3j;p -I- 

from the Ei term, where g{r) is a radial function to be chosen. This gives 

-(^mP) = E^gir) + qir) + 


2{r-2M)r 

( (r - 2M}q'(r) _ (r-2M)q"(r) 
V j ,2 2 r 


TT- 


(5.12) 
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By Lemma [HI the first terms are non-negative if and only if the inequalities 


0 < g{r) < min(^i/'(r) 


max 



Mfjr) 

(r — 2M)r 


_ Mf{r) f{r){r - 3M) \ 

(r —2M)r’ {r — 2M)r / 

- + y'ir)^ < q{r) < -g{r) 


(5.13a) 


Mfjr) 

(r — 2M)r 


+ y'{r), 
(5.13b) 


hold. We clearly see that /(r) must satisfy the condition f{r){r — 3M) > 0. This means that it 
must change sign at r = 3M. It is also important that it is bounded and vanishes at r = 2M, so 
we can dominate Jaa' with Haba'B'^^^ ■ A natural choice is therefore 


, (r-3M)(r-2M) 

/(■') =-- 


(5.14) 


Now, g{r) must be chosen such that 


n / . ^ ^ • ({r-3Mf 3M{r-2M)^ 

Q<g{r)<mni^ - ^-3 -,- — -j. (5.15) 

We do not want to saturate the inequalities except at r = 2M and r = 3M, and this can be 
achieved by setting 


g{r) = Cl 


3M(r-3M)2(r-2M) 
Ar^ ’ 


for some constant 0 < ci < 1. Then, the upper bound for q{r) is 

3M(r - 2M)(r2 - ci(r - 3Mf) 


qir) < 


4r® 


We can therefore set 


, , 9M^(r — 2M){2r — 3M) 

^^^ ^ 0. 


(5.16) 


(5.17) 


(5.18) 


A direct calculation shows that with /(r), g{r), and q{r) given by (15.141) . (I5.16L and (I5.18p . 
respectively, the inequalities (15.131) are satisfied everywhere and saturated only at r = 2M and 
r = 3M. We now choose ci = 5/6. With these definitions we get using (I5.12p . 


-(^i.iP) = Ai 


(r - 3M)2(14M2 - 7Mr + Ar^) M{90M^ - 105M2r -h 28Mr^ + r^) 
16r^ ’ 

(r — 3M)^(10M^ — 5Mr + Ar^) 


Aj-5 


16r® 


Paa' 


-h |%P)(r - 3M)2(r - 2M) 27M^{r - 5M){r - 2Mf 


8r5 


2r8 


W 


TT 


M\Pf\'^{r-3Mf{r-2M) ^ M|%|2(r - 2M)(r2-h 66Mr - QQM^) 
{\Px? + I/5yP)(^ - 3M)2(2r2 - 3Mr + 6M^) 

+ 4r5 

5M(|/3.f |2 |%P)(r - 3Mf{r - 2M) 27M^{r - 5M){r - 2Mf 


> g I/3 a A'In deg + 

27M^{r -5M){r -2Mf 

2/^8 


87-5 2r8 

5M{\l3f\^ + |/3gP)(r- - 3M)2(r - 2M) 


W 


TT 


8r5 


W 


TT- 


(5.19) 
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The last term in (|5.19l) is not positive everywhere and therefore must be estimated. Using the 
Hardy estimate in Lemma [TS] we get 



-{S)is'P)dnSr > 



2 

l,deg 


> 


M(r - 2M)(5r-3 - 8AMr^ + A28M^r - , 

+ ^^ 

/• 1 ,^ ,2 , ^ 

Jg + lOOr^ 

gl/^AA'Il.deg + 2 ^qq ^4 (5.20) 

□ 


6. Proof of the main theorem 

Proof of Theorem[^ First of all Lemma fTTl and the divergence theorem gives that E^{T, 2 ) = 
E^CSi). Lemma [TTl Theorem [T51 and the divergence theorem yield 

^ g l/5a|l,deg + ^qq ^4 l^af)|2'^Mf2 ^ Pad/lO = (Ui) — (U 2 ). (6-1) 

Corollary [TTl and the conservation of E^(E) then gives 

r ^ M 19 1 9 

8l/3al?.deg + Y^,\Zat\ldnn < - j^E,{E 2 ) = -Uj(Ei)- (6.2) 

□ 
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